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Expressions  are  derived  for  the  spectral  densities  $>T(K),  of  the 

temperature  and  velocity  structure  functions  of  atmospheric  turbulence, 
and  for  the  corresponding  Bom  approximation  far-field  acoustic  scattering 
cross-sections,  due  to  homogeneous  isotropic  stationary  ensembles  of 
self-similar  localized  turbules  having  many  different  scale  lengths.  It  is 
shown  that  for  some  range  i^min  <K<  Kmax/  the  "inertial  range,"  the 
spectral  densities  obey  power  laws  with  dependence  K~Pt,  K~Pv.  The 
exponents  (PT ,  Pv)  depend  only  on  choices  of  scaling  relations  and  are 
independent  of  turbule  morphology.  Only  Km-m,  Kmax,  and  the  values  of 
the  spectral  densities  outside  the  inertial  range  are  morphology- 
dependent.  Expressions  for  Kmin  and  Kmax  are  derived  in  terms  of  inner 
and  outer  scale  lengths  in  the  turbule  ensemble.  If  the  turbule  scale 
lengths  aQ  are  chosen  to  be  in  geometric  sequence  (aa+1/aa  =  constant 
independent  of  a),  and  if  the  power  law  is  given  as  PT  =  Pv  =  11/3,  the 
Kolmogorov  spectrum  in  the  inertial  range,  then  not  only  must  the  turbule 
velocity  and  temperature  amplitudes  scale  as  aj/3,  the  usual  result,  but 
also  the  turbule  packing  fractions  must  be  independent  of  scale  length. 
Expressions  for  the  structure  parameters  ((?£,  C%)  that  occur  in  the  usual 
Kolmogorov  spectra  are  obtained  in  terms  of  the  turbule  model 
parameters.  It  is  also  shown  that  quasi-Gaussian  spectra  result  for  the 
choice  PT  =  Pv  —  0  and  Gaussian  turbule  morphology. 
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1.  Introduction 


This  is  the  second  in  a  series  of  technical  reports  that  examine  a  turbule 
model  of  atmospheric  turbulence  and  the  acoustic  scattering  predicted  by 
the  model.  The  previous  report  [1]  laid  the  groundwork  for  this  and 
succeeding  reports  in  the  series.  In'  Goedecke  et  al  [1],  the  general 
properties  of  the  Bom  approximation  far-field  acoustic  scattering 
predicted  by  Monin's  equation  [2]  were  obtained.  In  addition,  expressions 
for  the  acoustic  scattering  amplitudes  and  cross-sections  were  derived  for 
individual  model  turbules  of  a  given  scale  length  and  for  their 
orientational  averages,  and  it  was  shown  that  an  ensemble  of  randomly 
oriented  turbules  of  arbitrary  morphology  may  be  replaced  by  an 
equivalent  ensemble  of  spherically  symmetric  nonuniformly  rotating 
turbules,  with  randomly  directed  rotation  axes. 

In  this  report,  we  connect  the  turbule  model  predictions  to  the  structure 
function  predictions  for  the  case  of  isotropic  homogeneous  fully 
developed  steady-state  atmospheric  turbulence.  This  is  an  essential  step 
toward  the  ultimate  goal  of  describing  acoustic  propagation  and  scattering 
in  inhomogeneous  anisotropic  turbulence  using  a  turbule  model.  In 
section  2,  we  construct  an  ensemble  of  self-similar  turbules  of  many 
different  scale  lengths,  each  with  random  location  in  a  bounded  volume  V 
and  with  random  orientation.  We  adopt  general  scaling  laws  in  which  the 
number  of  turbules  of  a  given  scale  length  aa  and  their  velocity  and 
temperature  variation  amplitudes  scale  according  to  powers  of  aa,  and  in 
which  the  spectrum  of  scale  lengths  follows  a  general  power  law  that 
includes  the  usual  fractal  geometric  sequence  in  which  aa+i/aa  is 
independent  of  aa.  We  also  show  that  a  kinetic  energy  cascade  model 
similar  to  that  of  Kolmogorov  yields  one  connection  among  the  scaling 
exponents.  We  develop  expressions  for  the  spectral  densities 
$t(K),  of  the  temperature  and  velocity  structure  functions  in 

terms  of  these  turbule  model  parameters.  We  show  that  a  power  law 
scattering  spectrum  generally  exists  for  some  range  Km[n  <K<  Kmax,  the 
"inertial  range,"  in  which  the  cross-sections  each  depend  on  K  to  some 
power  that  could  be  different  for  the  temperature  than  for  the  velocity 
scattering.  We  also  show  that  the  power  laws  are  independent  of  turbule 
morphology,  and  that  only  Km  m  and  Kmax  and  the  behavior  of  the  spectral 
densities  for  K  outside  the  inertial  range  are  morphology  dependent.  We 
derive  expressions  for  (Km jn,  Kmax)  in  terms  of  the  inner  and  outer  scale 
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length  of  the  ensemble,  and  the  characteristic  widths  of  the  spectral 
densities  of  individual  turbules.  We  show  that,  if  we  wish  to  obtain  a 
spectrum  in  which  both  spectra  depend  on  K~n^3  in  the  inertial  range, 
then  for  the  fractal  sequence  only,  not  only  must  the  velocity  and 
temperature  amplitudes  scale  as  alf3,  the  usual  result,  but  also  the  turbule 
packing  fractions  must  be  scale-invariant.  We  show  that  a  Gaussian 
scattering  spectrum  requires  quite  different  scaling  exponents  than  a 
r-11/3  spectrum.  We  also  obtain  expressions  in  terms  of  the  turbule 
ensemble  parameters  for  the  structure  parameters  (C|,  C$)  that  occur  in 
the  Kolmogorov  structure  functions.  We  investigate  the  specific  behavior 
of  the  spectral  densities  versus  K  for  two  example  turbule  structures. 

In  section  3,  we  express  the  Bom  approximation  far-field  cross-sections  as 
functions  of  scattering  angle  6  for  acoustic  scattering  by  the  velocity  and 
temperature  fluctuations  of  the  turbulence  derived  in  section  2,  and  show 
under  what  conditions  the  cross-sections  deviate  appreciably  from  a 
power  law  dependence  on  sin(0/2). 

Finally,  in  section  4  we  summarize  and  discuss  our  results  and  plans  for 
further  work. 


2.  Turbule  Model  of  Homogeneous  Isotropic  Steady-State 
Atmospheric  Turbulence 


2.1  Model  Turbules 

We  model  turbulence  contained  in  a  volume  Vs  as  an  ensemble  of 
self-similar  stationary  localized  turbules  of  different  scale  lengths.  On  the 
average,  in  Vs  we  allow  Na  turbules  of  scale  length  aar  a  =  (1,  Ns),  where 
Ns  is  the  total  number  of  different  scale  lengths  in  the  ensemble,  so  that 
N  =  ZaNa  is  the  total  ensemble  average  number  of  turbules  in  Vs,  in 
steady  state.  We  assume  in  this  report  that  the  aa  are  much  smaller  than 
the  scale  length  as  of  the  volume  Vs. 

We  take  ax  as  the  largest  scale  length  in  the  ensemble,  and  a^s  as  the 
smallest;  these  lengths  define  the  outer  and  inner  scales  of  the  turbulence, 
respectively.  Each  turbule  is  characterized  by  the  general  scalable  static 
temperature  variation  and  solenoidal  velocity  fields  (AT0(r),  v0(r))  chosen 
in  the  previous  report  in  this  series  [1].  For  turbule  number  n,  we  have 

von  (r)  =  V4n  x  An(0,  AT0n(r)  =  rn(£n)  (1) 

where 

=  (r  —  bn)/Oji)  (2) 

b„  is  the  "location"  of  the  turbule,  i.e.,  the  point  about  which  the  turbule  is 
localized,  and  An(£n)  is  an  unrestricted  vector  field. 

For  convenience,  we  allow  uncorrelated  locations  bn.  This  means  that  the 
joint  probability  distribution  for  the  locations  of  the  N  turbules  may  be 
written  as 

P(b, . b„)  =*,(■>,). ..j>„(b„),  (3) 

the  product  of  N  one-particle  distributions,  with 

J  d3bpn(b)  =  1,  alln.  (4) 

We  also  want  the  turbulence  to  be  homogeneous.  This  means  that  the 
ratios  of  the  number  densities  of  turbules  of  different  scale  lengths  should 
be  independent  of  position  in  Vs ;  this  requires 

Pn(b)=p(b),  (5) 
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independent  of  n.  If  we  also  want  the  turbulence  to  be  uniform  in  Vs,  then 
we  must  choose 

p(b)  =  V.-‘,b«V.  (6) 

=  0,  otherwise. 

An  isotropic  ensemble  consists  of  randomly  oriented  copies  of  turbules  of 
arbitrary  morphology  for  each  scale  length,  and  is  thus  characterized  by 
spherically  symmetric  envelope  functions  B2 ,  B2  as  defined  by  Goedecke 
et  al  [1],  namely. 


(|Tn(KaJ|2)  =  7t3(8Ta)2B2T(Kaa)  (7) 

(ini(Kaj)  =  0  (8) 

(AniA*nj(Kan))  =  1 7r3  Stj  v^B2{KaQ).  (9) 


The  Tn(yn)  and  Ani(yn)  are  the  Fourier  transforms  of  Tn(£n)  and  Ani (£„)> 
respectively,  as  defined  by  Goedecke  et  al  [1],  namely, 

Tn(y )  =  Jd^e-^TniS),  A  n(y)  =  Jd3^-^  An(£),  (10) 

where  the  integrals  are  extended  over  all  £-space,  since  the  individual 
turbules  are  localized  and  their  scale  lengths  an  are  assumed  to  be  much 
smaller  than  the  scale  length  as  of  Vs.  Here,  the  expectation  involves 
averaging  over  random  orientations,  so  the  amplitudes  ( STa ,  va)  and  the 
argument  of  the  envelope  functions  depend  only  on  the  scale  length  index 
a.  The  envelope  functions  themselves  are  independent  of  a;  that  is,  they 
are  the  same  functions  of  their  arguments  for  all  a.  This  provides 
self-similarity.  The  factors  of  it3  and  1/3  are  inserted  for  convenience  as  in 
Goedecke  et  al  [1]. 


2.2  Scaling 

In  order  to  describe  the  complete  ensemble,  we  assume  that  the  quantities 
(Na,  5Ta,va)  scale  with  aa.  We  put 
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where  (/?,  7,  v)  are  parameters.  In  addition,  we  must  decide  how  to  relate 
the  scale  lengths  to  the  index  a.  One  relation  that  has  been  used  is  [3] 


aa  =  ale  ^ Q  n  >  0,  (12) 

where  n  is  a  parameter  that  is  determined  by  Ns,  the  number  of  scale 
lengths,  and  the  ratio  (ojv3  /a  1 )  of  inner  to  outer  scale  lengths: 

H  =  —(Ns  —  l)-1  lnm,  m  =  a^Ja1.  (13) 

Equation  (12)  implies  that  the  scale  lengths  form  a  geometric  sequence,  in 
which 

^2/^1  =  ^3/^2  =  •  •  •  =  ®  ^*5 

that  is,  the  ratio  of  successively  smaller  scale  lengths  is  a  constant  whose 
value  lies  between  zero  and  unity.  This  is  a  kind  of  fractal  scaling  [3]. 

A  general  power  law  scaling  relation  is  given  by 

aa  =  (1  +  qfJi(a  -  q>  0,  fi  >  0,  (14) 

where  here  [i  is  determined  in  terms  of  (m,  q)  by 

M  =  «-1(m-«-l)/(^-l).  (15) 

Equations  (14)  and  (15)  actually  reduce  to  equations  (12)  and  (13) 
respectively,  when  q  goes  to  0;  so  in  what  follows  we  may  use  just 
equations  (14)  and  (15),  with  q  >  0.  We  will  also  make  use  of  the 
Kolmogorov  concept  of  energy  transfer  in  fully  developed  (steady-state) 
turbulence,  in  which  energy  is  provided  to  the  ensemble  at  the  largest  or 
outer  scale.  The  largest  eddies  continually  form  but  are  unstable  because 
of  their  large  Reynolds  numbers,  and  thereby  continually  fragment  into 
smaller  eddies,  which  in  turn  fragment  further,  etc.  This  cascade  continues 
down  to  eddies  of  a  size  small  enough  to  be  stable,  that  is,  to  eddies  whose 
Reynolds  numbers  are  of  order  unity.  These  smallest  or  inner  scale  eddies 
dissipate  almost  all  the  energy  that  is  being  input  at  the  largest  scale.  In 
steady  state,  all  the  ensemble  average  quantities  (Na,  va,  6TQ),  and  the 
relationships  like  equations  (12)  or  (14),  are  constant  in  time.  By 
dimensional  analysis  and  from  the  fluid  equations,  the  kinetic  energy 

transfer  rate  £Ka  from  turbules  of  scale  length  aQ  to  the  next  smaller  is 
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£Ka=  0 C)(Na)(al)(6Ta)(vl/aQ ), 


(16) 


where  C  is  a  constant  with  dimension  of  reciprocal  volume.  That  is,  the 
energy  transfer  rate  is  proportional  to  the  number  of  turbules  of  size  aa  in 
V ,  the  volume  of  each,  the  kinetic  energy  per  unit  mass  v^/2  of  each,  and 
the  characteristic  rate  of  transfer  va/aa.  The  Kolmogorov  model  consists  of 
neglecting  dissipation  in  all  eddies  except  the  smallest.  This  means  that,  in 

steady  state,  £Ka  is  independent  of  a,  for  a  =  (1,  Ns-1).  Using  equation 
(11)  in  this  model  then  yields  one  ratio  among  the  parameters  ((3,  v)\ 

£  =  3^  +  2.  (17) 

In  the  atmosphere,  the  ratios  lAT^/T^  and  v^jcx  are  usually  of  the  same 
order.  This  implies  that  our  turbule  temperature  variation  amplitudes  STa 
should  be  proportional  to  va,  whereby  we  get  directly  from  equation  (11) 

7  =  i/.  (18) 

On  the  other  hand,  if  we  apply  the  Kolmogorov  energy  cascade  model  to 
the  thermal  energy  constant  of  the  turbules  in  our  ensemble,  we  get  for  the 
thermal  energy  transfer  rate  £Ta  from  turbules  of  one  scale  length  to  the 
next  smaller 


£Ta=  (C')(Na)(al)(STa)(va/aa),  (19) 

where  C'  is  a  constant.  This  says  that  the  energy  transfer  rate  is 
proportional  to  the  thermal  energy  content  of  a  turbule,  which  is 
proportional  to  STa.  If  we  require  this  rate  to  be  independent  of  a,  then  we 
get  from  equations  (11)  and  (14) 


7  =  2i/.  (20) 

This  is  equivalent  to  stating  that  6Ta  is  proportional  to  v\.  We  see  in  the 
next  section  that  these  two  possibilities  of  equations  (18)  and  (20)  yield 
quite  different  behavior  for  the  velocity  and  temperature  variation 
structure  functions  of  the  ensemble. 

2.3  Structure  Function  Spectra 
2.3.1  General 

The  velocity  and  temperature  variation  structure  functions  are  important 
,  quantities  that  characterize  atmospheric  turbulence  [4,5].  In  this  section, 


we  obtain  formulas  for  the  spectral  densities  of  these  structure  functions 
from  our  turbule  model,  for  isotropic  homogeneous  turbulence. 

The  temperature  variation  and  velocity  fields  of  the  turbulence  in  Vs  are 
just  the  superposition  of  those  of  the  localized  turbules  in  Vs: 

AT.M  =  EAI^r),  =  (21) 

n  n 

We  are  interested  in  the  Fourier  transforms 

Af0(K)=  f  d3re~iK  r AT0(r),  vq(K)  =  f  d3re~iK-rv0(r).  (22) 

Specifically,  we  wish  to  obtain  general  expressions  for  the  spectral 
densities 

*Tm  =  (|Ar„(K)|2),  *v(K)  =  (s„(K)s;(K)),  <23) 

where  the  expectations  ()  imply  averaging  over  random  turbule  locations 
b  and  random  turbule  orientations,  as  discussed  in  section  2.1  above. 

Combination  of  equations  (1),  (2),  (10),  (21),  and  (22)  yields 

^(K)  =  *X!e",K'b’,OnKxAtl(Kan)  (24) 

n 

Af0(K)  =  (25) 

n 

The  spectral  densities  of  equation  (23)  are  then 
$t(K)  =  £a$(|f„(Ka„)|2} 

n 

+  lp(K)|2Efln(r„(Kan))^a3(f;(Kan))  (26) 

"  t^n 

—  ^ipqtjrsKpKr  ^  an  ^-4„9(Ka„)A*  ,.(Kan)^  .  (27) 

n 

where  equation  (8)  was  used.  Here  eipq  is  the  Levi-Civita  symbol,  and 

p(K)  =  [  dzbe~iK  hp{ b)  (28) 

•'Vs 
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is  the  Fourier  transform  of  the  "one-particle"  location  distribution  of 
equation  (6). 

In  the  standard  treatment  [4],  the  mean  temperature  in  Vs  is  assumed 
equal  to  the  remote  reference  temperature  T^,  or  /v/rAT0(r)  =  O.We 
will  assume  that  here,  in  order  to  facilitate  comparisons  of  our  turbule 
model  spectra  with  the  commonly  used  spectra.  For  simplicity  we  assume 
that  in  the  ensemble, 

(fn(Kan))  =  ±irV26TaBT(Kaa)  =  ±  (|  fn(Kan)|2)^  ,  (29) 

with  equal  numbers  having  the  (+,  — )  signs.  This  ensures  that 
f  d3rAT0  (r)  =  0.  It  also  takes  advantage  of  the  result  of  Goedecke  et  al 
[1],  that  an  isotropic  ensemble  of  a  given  scale  length  may  often  be 
replaced  by  an  ensemble  of  spherically  symmetric  turbules.  For  such 
turbules  with  a  positive  definite  envelope  function,  equation  (29) 
automatically  would  be  valid.  With  this  assumption,  we  have 

£>3  (rn(Kan))  =  0.  (30) 

n 

Then,  from  equations  (7),  (9),  and  (26)  to  (30),  we  get 

$r(K)  =  (l  —  |p(K)|2)  7r3^iVQ((5Ta)2a®S2(KaQ)  (31) 

n 

$V.(K)  =  (Sij  —  KiKj)(ir3/3)^2Nav2a^(Kaa)2B2(Kaa).  (32) 


Note  that  the  sums  over  turbules  in  equations  (26)  and  (27)  are  replaced  by 
sums  over  size  index  a  in  equations  (31)  and  (32),  with  Na  included. 

The  factor  involving  |  p(K)|2  in  equation  (31)  comes  from  writing 
En  E^n  =  En  E^  -  En  in the  second  term  of  equation  (26),  and  using 
equations  (29)  and  (30).  In  general,  if  fy  d3rAT0  (r)  =  0,  then  $r(K  =  0) 
must  be  zero;  the  factor  (1  —  |  p(K)  |2)  in  equation  (31)  ensures  this,  since 
p( K  =  0)  =  1,  because  p(b)  is  a  probability  density  (see  eq.  (28)). 

It  is  important  to  note  that  p(K)  is  extremely  small  except  for  K  near  zero. 
For  example,  suppose  that  14  is  a  spherical  volume  of  radius  as  centered  at 
the  origin.  Then,  from  equations  (6)  and  (28) 

p(K)  =  3 (Kas)~3  [sin  Kas  —  Kas  cos  Kas } .  (33) 
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Thus  for  large  Kas,  | p(K)|2  is  smaller  than  (Kas)~4  «  1.  Therefore,  in 
the  spectral  density  $T(K),  the  |  p(K)|2  factor  may  be  dropped  out,  except 
for  very  small  K  such  that  Kas  <  1.  But  its  presence  is  essential  in  order  to 
comply  with  the  assumption  that  fyd3rAT0(r)  =  0. 

This  assumption  is  not  necessary;  if  it  is  not  made,  then  the  spectrum 
$t(K)  will  have  a  different  behavior  for  K  -*  0  than  that  of  equation  (31). 
We  will  investigate  this  in  a  later  report. 

Note  that  the  velocity  spectrum  $?-(K)  includes  the  factor  of  (5{j  —  KiKj), 
characteristic  for  solenoidal  velocities  [4].  Also  note  that,  for  the  bounded 
envelope  functions  B2  that  will  be  used,  the  factor  ( Kaa )2  ensures  that 
$ij(K  — >  0)  — *•  0,  which  is  also  a  requirement  for  solenoidal  velocities, 
directly  related  to  the  results  of  Goedecke  et  al  [1]  and  Monin  [2]  that 
forward  acoustic  scattering  due  to  such  turbulent  velocities  is  zero. 

2.3.2  Inclusion  of  Scaling 

In  what  follows,  we  replace  the  sum  over  size  index  a  by  an  integral.  This 
will  be  valid  if  the  number  of  different  scale  lengths  in  the  ensemble  is 
large,  as  we  shall  assume.  We  then  write 

r Ns  raNs 

Y,  =  /  da  =  da /(da/ da)  \  da/da  =  fj,a~qa1+q  (34) 

a=l  •'i  ■'“i 

where  the  last  equality  results  from  equation  (14). 

From  equations  (31),  (32),  and  (34)  and  the  scaling  relations  (11)  to  (15),  we 
then  get  the  following  expressions  for  the  spectral  densities: 

&t(K)  =  (1  —  |p(K)|2)  (ir3N1(ST1)2a^/fj,)x~Pr^_i(mx,x)  (35) 

$V.(K)  =  (5ij  -  yky^N.vyj^x-^r^  (mx,  x),  (36) 

where  (m,  n)  are  defined  by  equations  (13)  to  (15),  and 

x  =  Ka1 ;  (37) 

Jj'v(mx,x)=  f  dyysB2Tv(y)  (38) 

Jmx 

PT  =  6  +  27  -  /?  -  q,  Pv  =  6  +  2i/  -  /?  —  q.  (39) 

Here  q  >  0;  q  =  0  corresponds  to  the  fractal  scaling  of  equation  (11). 
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Thus  the  spectra  of  equations  (35)  and  (36)  depend  on  the  integrals  of  the 
generic  form  of  equation  (38).  It  is  important  to  determine  qualitatively 
how  these  integrals  depend  on  (s,  m,  x).  In  all  cases,  the  envelope 
functions  B2  v  ( y )  go  to  zero  rapidly  for  large  y;  otherwise,  individual 
turbules  would  not  be  localized.  Also,  we  expect  the  B2(y)  to  be  bounded 
everywhere.  Thus,  if  x  is  large,  then  Js(mx,  x)  «  Js(mx,  oo)  to  very  good 
approximation,  as  long  as  the  B2(y)  go  to  zero  faster  than  y~s~l  for  large 
y.  Similarly,  if  mx  «  1,  then  Js(mx,  x )  «  Js( 0,  x)  to  very  good 
approximation,  as  long  as  the  ysB2(y)  — *•  yl,  with  t  >  —1,  for  y  — »  0.  Thus, 
in  many  cases,  there  will  exist  a  range  of  x  such  that 

Js(mx,x)  ~  Js(0,  oo)  =  constant  (40) 

to  a  very  good  approximation. 

For  this  range  of  x,  equations  (35)  and  (36)  show  that  the  spectra  have  a 
power  law  dependence  on  x  and  thus  on  K,  with  powers  (— PT,  —Pv). 
Conventional  language  defines  this  range  of  x  or  K  as  the  "inertial  range." 

We  note  that,  if  the  powers  (PT ,  Pv)  are  to  be  the  same,  then  from  equation 
(39),  we  must  have  equation  (18), 


7  =  i/,  (41) 

which  is  what  resulted  from  the  discussion  preceding  equation  (18),  not 
from  that  preceding  equation  (21).  We  shall  adopt  7  =  v.  Then  equations 
(17)  and  (39)  yield 

P  =  PT  =  Pv  =  4-v-q.  (42) 


2.3.3  Inertial  Range  Boundaries 

The  boundaries  of  the  inertial  range  of  x  may  be  estimated  as  follows.  The 
integrand  Is(y)  of  Js(mx,  x)  of  equation  (38)  is  y2B2(y );  for  s  >  0  and 
B2{y)  that  decrease  monotonically  faster  than  y~s  as  y  increases,  it  has  a 
single  maximum  at  y  =  ysm  given  by  /'(ysm )  =  0,  or 

Z/sm  =  3 sB(ysm)/ B  ( ysm )•  (43) 

The  integrand  has  a  characteristic  width  that  also  depends  on  s  and  the 
envelope  function  B(y).  We  define  ys±  by  requiring  that  Is(ys± )  be  some 
fraction  of  I(ysm)-  In  this  report  we  choose 

I^s(j/s±)  =  6  Is {ysm)>  Vs—  Vsm  ^  Vs+  ■  (44) 
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Then  it  is  clear  that  Js(mx,  x )  «  Js( 0,  oo)  for  mx  <  ys-  and  x  >  ys+.  That 
is,  essentially  pure  power  law  spectra  obtain  for  values  of  x  that  lie 
between  xs  min  and  xs  max,  where 

2-s  min  ~  2/s+>  Xs  max  ~  2/s— /m-  (45) 

If  x  <  xs  mm  or  x  >  xs  max,  the  dependence  of  the  spectra  departs 
significantly  from  that  of  a  power  law.  Note  (armin,  xmax)  for  $T(K)  are  in 
general  different  than  for  $?•,  because  s  is  different  for  the  two  cases  for 
PT  =  Pv  ,  and/or  because  the  B{y)  may  be  different.  Also  note  that  if 
xa  min  >  xs  max  in  equation  (45),  then  there  is  no  inertial  range. 

Morphology  Dependence 

It  is  important  to  examine  the  effects  of  changing  the  envelope  functions 
B2{y)  that  appear  in  the  integrals  Js(mx,  x)  of  equation  (38). 

We  may  change  an  envelope  function  in  several  ways.  One  way  is  merely 
to  alter  its  amplitude.  But  that  is  trivial,  because  the  previous  values  of  the 
functions  $T(K),  §\j(K)  could  be  maintained  by  commensurately  altering 
the  unknown  amplitudes  in  equations  (35)  and  (36). 

Another  way  is  to  replace  B(y)  by  the  same  functions  of  a  "stretched" 
argument,  B(y)  — >  B(ay).  This  is  equivalent  to  changing  the  spectrum  of 
scale  lengths  in  the  ensemble,  such  that  aQ  —ya'Q  =  aaQ,  and  altering 
(<5Tj ,  Vj )  appropriately  to  keep  the  same  values  of  $T(K),  ( K )  in  the 

inertial  range.  But  the  boundaries  of  the  inertial  range  will  be  shifted;  that 
is,  from  equation  (45),  we  get 

■^pmin  =  yp+/act  =  Kpmm/ot,  ^pmax  =  Bpmax/a.  (46) 
Outside  the  inertial  range,  the  spectra  will  be  changed. 

Another  way  is  to  replace  a  chosen  B(y)  by  a  different  functional  form 
B'(y).  Clearly,  this  will  change  the  boundaries  of  the  inertial  ranges,  in 
general,  but  inside  the  inertial  ranges,  (STt ,  v1 )  can  be  altered  to  preserve 
the  previous  values  of  the  spectra. 

Therefore,  we  may  conclude  the  following:  the  power  law  spectra  in  the 
inertial  ranges  are  completely  insensitive  to  all  changes  in  turbule 
morphology,  that  is,  alterations  of  the  envelope  functions  B(y).  Changes 
in  the  spectrum  of  scale  lengths  and/or  in  the  functional  form 
(morphology)  of  the  B(y)  irreducibly  influence  only  the  boundaries 

(-Kprmn,  Kp 

max  )  of  the  inertial  range  and  the  behavior  of  $T(K), 
outside  the  inertial  range. 

We  will  provide  graphs  of  numerical  examples  in  section  2.4.3. 
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2.4.1  Kolmogorov  Spectra 

From  experimental  and  scaling  considerations  [4],  the  power  law  part  of 
the  spectrum  is  expected  to  go  like  x~n^3r  so  P  =  11/3.  If  this  dependence 
were  valid  for  all  x,  the  spectral  densities  would  be  said  to  exhibit  the 
"Kolmogorov  spectrum." 

From  equations  (42)  and  (17),  this  requires 


u  =  l/3-q,  (3  =  3(1  -  q).  (47) 

We  know  that  q  >  0,  since  we  insisted  that  aa  get  smaller  as  a  increases. 
Also,  we  expect  that  Na  increases  as  the  scale  length  decreases;  from 
equations  (11)  and  (47),  this  requires  q  <  1.  It  is  especially  interesting  to 
note  that,  for  the  fractal  scaling  (q  =  0),  which  is  often  assumed  [2], 
equation  (47)  yields 


v  =  1/3  /?  =  3.  (48) 

The  v  =  1/3  result  follows  from  the  usual  energy  cascade  model,  upon 
requirement  that  the  kinetic  energy  transfer  rate  per  unit  mass  be 
independent  of  scale  length  [4].  The  usual  cascade  model  does  not 
consider  the  number  of  eddies  of  each  scale  length,  in  contrast  to  our 
model  of  section  2.2.  It  is  remarkable  that  only  with  fractal  scaling  does 
our  model  yield  not  only  v  =  1/3,  but  also  /3  =  3,  which  corresponds  to 
turbule  packing  fractions  Naa^/Vs  independent  of  scale  length. 

The  standard  structure  function  description  of  isotropic  homogeneous 
fully  developed  turbulence  finds  by  dimensional  analysis  that  the 
temperature  and  velocity  structure  functions  of  the  turbulence  must  be 
given  by 


DT(r)  =  ^(ro(rJ-T0(r2))^ 

)  =  cy/\ 

(49) 

DVTT{r)  =  rifj((v0i{ rj  -  v0t(r2))  (%(rj  ■ 

-»«(*,))) =cy/3, 

(50) 

where 

III 

1 

to 

II 

jjT 

r  =  r/r, 

(51) 

and  (C3,C%)  are  the  so-called  structure  parameters.  These  are  valid  for 
some  range  of  r  called  the  "inertial  range"  between  the  inner  and  outer 
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scale  lengths  (ajv,,  oj.  Assuming  (incorrectly,  of  course)  that  these  forms 
are  valid  for  all  r,  the  following  spectral  densities  are  easily  derived: 


$t(K)  =  8.19  C*VsK~ntz  (52) 

%{K)  =  (Sij  -  KiKj)(15.0  C2vVsK~n/z).  (53) 

It  is  clear  that  these  are  incorrect  for  K  — >  0,  but  they  are  valid  for  most  K. 
Comparing  these  forms  with  equations  (35)  and  (36),  putting  PT  =  Pv  = 
11/3,  neglecting  |  p{ K)  j2  <<  1  for  K  not  near  0,  and  taking  the  values  of 
the  integrals  of  equation  (33),  in  the  inertial  range  of  K,  we  get  the 
following  expressions  for  the  structure  parameters: 


C2t  =  (3.78/ p)(N1a^/Vs)(ST1  /a\^3)2  Jq/3(0,  oo),  (54) 

Cl  =  (0.69/ y)(Nla\/Vs)(v1  /alJz)2  oo).  (55) 

Thus  we  have  connected  the  structure  parameters  that  are  believed  to 
characterize  the  Kolmogorov  spectrum  of  isotropic  homogeneous 
turbulence  to  the  parameters  of  our  scaled  self-similar  turbule  model,  an 
essential  step.  Note  that  for  fractal  scaling  (q  =  0),  the  factor  (N^/Vg) 
may  be  replaced  by  (Naa^/Vs)  for  any  scale  length  a,  and  so  may  the 
factors  (5Tl/alJz,  Vj/aJ/3).  This  is  not  the  case  for  q  ^  0. 

2.4.2  Quasi-Gaussian  Spectrum 

Gaussian  spectra  involving  a  single  scale  length  have  been  used  fairly 
often  [5,6];  we  illustrate  this  for  an  isotropic  ensemble  of  turbules  of  a 
given  scale  length  in  Goedecke  et  al  [1].  It  is  interesting  to  note  that  our 
turbule  ensemble  model,  containing  many  scale  lengths,  allows  spectra  for 
$t(K),  that  are  Gaussian  for  most  K,  for  special  choices  of  the 

envelope  function  In  particular,  if  in  equations  (35)  and  (36)  we 

put 


PT  =  Pv  =  0,  B2T(y)  =  y2e-y^ 2,  B2v(y)  =  e^2,  (56) 

we  get 

$t{K)  =  (1  -  |p(JT)|2)(constant)(e-m212/2  -  e"*2),  (57) 

=  (5ij  -KiKj) (constant) (e-™2*2/2  _  e"*2).  (58) 
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For  K  — *•  0  and  for  K  — *•  oo,  these  go  to  zero,  as  they  should;  but,  for 
x  —  Kal  >>  1,  they  go  to  simple  Gaussians  involving  the  inner  scale 
length  clns,  since  mx  =  Ka^s.  These  are  a  special  case  of  a  power  law 
spectrum  with  P  =  0  for  the  "inertial  range"  of  K,  that  is,  for  mx  «  1 
and  x  »  1. 

From  equations  (17)  and  (39),  the  scaling  exponents  must  then  satisfy 

0  =  14  —  3 q  >  0,  '  v  =  7  =  4  —  q.  (59) 

If  we  were  to  choose  v  =  1/3  as  in  the  Kolmogorov  spectrum  with  fractal 
scaling,  this  again  yields  0  =  3,  but  requires  q  =  11/3,  a  rather  steep  power 
law  for  the  length  scales  in  the  ensemble  (see  equation  (14)). 


2.4.3  Examples 

It  is  important  to  give  some  examples  using  specific  envelope  functions,  in 
order  to  illustrate  some  of  the  results  of  the  last  several  sections.  We  will 
do  this  for  the  (quasi)  Kolmogorov  spectra,  having  dependence  x~n/3  in 
the  inertial  range. 

In  Goedecke  et  al  [1],  we  considered  two  example  envelope  functions,  a 
Gaussian  and  the  Fourier  transform  of  an  exponential: 

B2(y )  =  e-y*'2,  B2(y )  =  (1  +  y2/12)~6.  (60) 

The  factor  12  in  B2(y)  ensures  that  turbules  of  the  same  scale  length  have 
the  same  RMS  radius  [1].  For  P  =  11/3,  we  consider  the  normalized 
spectral  densities 


FT(x)  =  x  n/3J8/3(mx,x)/J8/3(0,oo),  (61) 

Fv(x)  =  x~n/3Ji4/3(mx,  x)/ Ji4/3(0,  oo),  (62) 

in  which  we  use  the  same  B(y)  for  both  temperature  and  velocity  spectra, 
for  convenience.  These  functions  are  the  factors  in  equations  (35)  and  (36) 
that  determine  the  boundaries  of  the  inertial  ranges  of  K  and  the  behavior 
of  $t(K),  outside  the  inertial  ranges. 

From  equation  (38),  we  have 

Jjl(mx,x)  =  f  dyyse~y2^2,  (63) 

Jmx 

J® (mx,  x)  =  f  dyys(l  +  y2/12)~6.  (64) 

Jmx 
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These  integrals  were  evaluated  for  s  =  (8/3, 14/3),  for  m  =  (10~3, 10-4), 
which  are  realistic  values  for  the  ratio  aN/ax  of  inner  to  outer  scale  length. 
That  is,  inner  scale  lengths  may  be  of  the  order  of  millimeters,  while  outer 
scale  lengths  may  be  of  the  order  of  tens  to  hundreds  of  meters. 

Figures  1  and  2  are  plots  of  log  F£  and  log  versus  log  a:  for  m  =  10-3 

and  10"4,  respectively;  figures  3  and  4  are  plots  of  log  F$  and  log  F Jf 


Figure  1.  Normalized 
isotropic  homogeneous 
ensemble  temperature 
spectra  for  m  =  10-3  as 
functions  of  outer  scale 
size  parameter  ka1 . 


Figure  2.  Normalized 
isotropic  homogeneous 
ensemble  temperature 
spectra  for  m  —  1CT4  as 
functions  of  outer  scale 
size  parameter  kax . 
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versus  log  x  for  m  =  10“3  and  10-4,  respectively.  The  plots,  of  course, 
coincide  in  the  common  portions  of  their  inertial  ranges;  as  discussed 
earlier,  this  coincidence  can  always  be  achieved  for  the  actual  spectra,  for 
any  choices  of  the  envelope  functions,  by  adjusting  the  (unknown) 
parameters  (STX ,  v1 ).  But  the  exponential  and  Gaussian  envelopes  yield 
slightly  different  inertial  range  boundaries,  quite  different  behavior  for 
x  >  xs  max/  and  the  same  behavior  but  different  values  for  x  <  xsm[n. 


Figure  3.  Normalized 
isotropic  homogeneous 
ensemble  velocity  spectra 
for  m  =  10“3  as  functions 
of  outer  scale  size 
parameter  kax . 


10°  102  104 


Outer  scale  size  parameter 


Figure  4.  Normalized 
isotropic  homogeneous 
ensemble  velocity  spectra 
for  m  =  1CT4  as  functions 
of  outer  scale  size 
parameter  kax . 
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We  may  use  the  results  of  section  2.3.3  to  estimate  the  inertial  range 
boundaries,  and  compare  the  estimates  with  figures  1  to  4.  Using 
equations  (43)  and  (60),  we  get 

Vim  =  sV2>  yesm  =  «1/2/(l  -  «/12)1/2  (65) 

for  the  (Gaussian,  exponential)  envelopes,  respectively.  Then  by  numerical 
-solution  of  equation  (44),  using  equation  (45),  we  obtain  the  following 
inertial  range  boundaries  for  the  temperature  spectra,  with  s  =  8/3: 


^max  ~  0.49/m  (66) 

^max  «  0.52/m.  (67) 

For  the  velocity  spectra,  with  s  =  14/3,  the  corresponding  boundaries  are 


xz 


!  3.19, 
:  4.45, 


^ax  ~  0.94/m  (68) 

Xmax  ~  1-08/m.  (69) 

These  results  compare  reasonably  well  with  the  boundaries  apparent  in 
figures  1  to  4. 

Based  on  the  discussion  in  section  2.4,  if  we  stretch  the  argument  of  the 
envelope  functions  {B(y)  — »  B(ay ))  or,  equivalently,  change  the  length 
scales  (aQ  — +  aaa),  then  the  boundaries  (i^min,  Kmax)  of  the  inertial  ranges 
of  K  change  according  to  equation  (46).  That  is,  if  we  knew  or  assumed  a 
value  of  ax,  then  we  would  get  from  x  =  Ka1  and  equation  (66) 

^min  ~  3-19/Oj,  w  0.49/mOj  =  0.49/ma^  (70) 

and  similarly  from  equations  (67)  and  (69).  If  we  then  change  a1  by 
<2j  — *  aaj,weget 


.9 

min 

,e 

min  ' 


:  3.70, 
:  6.34, 


•^min  ^  3.19/aOj  =  K^Ja,  etc. 


(71) 
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3.  Acoustic  Scattering 


For  an  isotropic  homogeneous  turbulence  in  a  volume  VT,  the  far-field 
Bom  approximation  differential  cross-section  for  scattering  of  an  acoustic 
plane  wave  is  given  by  [4] 

a(r)  =  aT  (? )  +  av  (? ) ,  (72) 

where 

aT  (r)  =  (/c2/47rT00)2cos2  6$T(K)  (73) 

av(  f)  =  (k2 /2nCoo)2  cos2  OkikjQVjiK)-,  (74) 

and  here, 

K  =  fcf  -  k,  K  =  |K|  =  2fcsin(0/2),  (75) 

where  r  is  the  observation  direction,  k  is  the  incident  plane  wave 
propagation  vector,  with  k  =  2n/\,  A  =  wavelength,  6  is  the  scattering 
angle,  the  angle  between  r  and  k,  0  <  0  <  7 r,  and  $T(K)  and  ^(K)  are 
the  spectra  of  the  turbulent  temperature  variations  and  velocity, 
respectively,  given  by  equations  (35)  and  (36).  In  equation  (73),  Too  is  the 
reference  temperature  outside  the  volume  VT,  which  has  been  chosen 
equal  to  the  mean  temperature  inside  VT,  as  discussed  in  the  previous 
section.  In  equation  (74),  we  have 

kikj^iji K)  =  kikj  ( Sij  -  kikj )  GV(K)  =  cos2{6 /2)GV{K),  (76) 

where,  from  equation  (36), 

GV(K )  =  (tt3^  vyj^x-P'’  JvPv+1  (mi,  x),  (77) 

and  from  equations  (37)  and  (75) 

x  =  Kax  =  2 kar  sin(0/2).  (78) 

Equations  (73)  and  (74)  are  valid  for  cases  in  which  VT  is  the  scattering 
volume  Vs,  or  in  which  Vs  <VT,  but  the  scale  length  as  of  Vs  satisfies 
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as»  aQ,  as  discussed  by  Goedecke  et  al  [1],  and  in  equations  (35)  and 
(36),  JVj  is  the  mean  number  of  turbules  of  scale  length  in  Fs. 

The  cross-sections  in  equations  (73)  and  (74)  are  thus  proportional  to  the 
normalized  spectral  densities  in  equations  (61)  and  (62),  and  figures  1  to  4 
therefore  reveal  the  crucial  behavior  of  the  scattering  cross-sections 
aT  (?),  av(r)  for  the  sample  envelope  function  considered  in  section  2.4.3, 
except  for  the  factors  of  cos2  9  and  cos2  6  cos2(0/2).  It  is  important  to  note 
that  here,  the  boundaries  of  the  inertial  range  of  x  or  K  translate  into 
boundaries  of  the  inertial  range  of  scattering  angles  (0min,  #max)-  For 
example,  for  the  Gaussian  envelope  function  and  the  x-11/3  power  law  in 
the  inertial  range,  we  have  from  equation  (66)  the  temperature  scattering 

sin(0min/2)  «  (3.19)/(2 fcaj,  sin(0max/2)  «  (0.49)(2mfca1 ).  (79) 

Clearly,  sin(0/2)  <  1  for  0  <  6  <  ir.  So,  if  mkal  is  less  than  or 
approximately  equal  to  0.25,  the  upper  inertial  range  boundary  is  never 
reached;  that  is,  the  power  law  spectrum  obtains  for  scattering  angles  out 

tO  6  =  7T. 

For  example,  suppose  the  following  reasonable  values  are  chosen: 
m  =  10-4,  a1  =  102m,  A  =  0.314m  — ♦  k  «  20m-1,  ka  «  2  x  103.  Then,  from 
equation  (79) 


sin(0min/2)  «  1.6  x  10  3,  sin(0max/2)  w  2.5.  (80) 

So  the  power  law  spectrum  would  obtain  for  almost  all  angles;  the  effect 
of  the  factor  1  —  |p(K)|2  in  <&T(K)  might  dominate  at  angles  as  small  as 
this  #min-  If  ka1  were  much  smaller,  then  0m;n  would  be  observable,  but  not 
9max)  if  kal  were  larger,  then  0max  would  be  observable  in  principle,  but 
0mm  would  be  too  small  to  be  observed. 
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4.  Summary  and  Discussion 


In  this  report,  we  have  established  that  a  self-similar  ensemble  of  localized 
stationary  turbules  of  arbitrary  morphology  and  many  different  scale 
lengths,  but  with  random  orientations  and  locations  within  a  volume  14, 
predict  spectral  densities  (<f>T(K),  $?•(#))  of  the  temperature  and  velocity 
structure  functions  of  the  turbulence  that  have  the  following  properties: 

1.  The  spectra  obey  power  laws  {K~Pt,  K~Pv)  in  ranges 

Km in  <K<  Kmax,  conventionally  called  the  inertial  ranges,  which 
may  be  different  for  4>T(K)  than  for  ^(K).  The  values  (PT,PV)  are 
determined  by  choices  of  scaling  exponents,  and  are  independent  of 
turble  morphology,  as  are  the  values  of  the  spectral  densities  in  the 
inertial  range.  The  boundaries  (Kmin,  -Kmax)  of  the  inertial  range  and 
the  behavior  of  the  spectral  densities  for  K  outside  the  inertial  range 
are  sensitive  to  turbule  morphology. 

2.  The  choice  PT  =  Pv  =  11/3,  corresponding  to  the  Kolmogorov 
spectrum  with  a  fractal  size  scaling  aa+i/aQ  =  constant,  yields  the 
usual  scaling  relation  vaal/ 3  =  constant,  and  also  requires  that 
packing  fractions  of  turbules  in  14  are  independent  of  scale  length. 
Expressions  for  the  structure  parameters  (C^,  C%)  of  the  Kolmogorov 
spectra  were  obtained  in  terms  of  turbule  model  parameters. 

3.  The  choice  Pv  =  PT  =  0  yields  a  quasi-Gaussian  spectrum  for  a 
Gaussian  ensemble  average  turbule  morphology. 

4.  The  first  Bom  approximation  far-field  acoustic  scattering 
cross-sections  due  to  the  turbulent  temperature  and  velocity 
fluctuations  exhibit  a  power  law  dependence  (sin#/2)_Pr, 
(sin0/2)_Pv  in  inertial  ranges  6m jn  <  0  <  0max  determined  by 
(Kmin,  Kmax),  but  deviate  markedly  from  this  dependence  for 
scattering  angles  6  outside  the  inertial  ranges.  Depending  primarily 
on  the  values  of  the  acoustic  wavelength  and  the  outer  and  inner 
scale  lengths  of  the  turbulence,  these  deviations  may  not  be 
observable  in  practice,  because  the  power  law  may  be  valid  from 
very  small  angles  out  to  6  =  180°. 

There  is  quite  a  bit  more  that  can  be  done  with  turbule  models  of 
atmospheric  turbulence.  For  example,  treatment  of  cases  in  which  the 
mean  temperature  in  14  is  not  equal  to  the  background  reference 
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temperature  can  be  done;  this  will  yield  a  different  power  law  for  $T(K) 
for  small  K  than  for  large  K.  Also,  consideration  should  be  given  to 
situations  that  must  often  occur  in  practice,  in  which  the  scattering  volume 
Vs  is  smaller  than  some  of  the  large-scale  turbules,  and/or  the  observation 
distance  is  not  in  the  far  field  of  Vs.  These  situations  may  change  the 
results  significantly.  They  will  also  be  considered  in  future  reports. 

The  ultimate  goal  of  the  research  reported  in  this  series  of  reports  is  to 
describe  acoustic  scattering  and  propagation  in  anisotropic  and/or 
inhomogeneous  turbulence.  It  is  hoped  that  a  turbule  approach  will  be 
particularly  appropriate  for  this  description. 
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